Positive P simulations of spin squeezing in a two-component Bose condensate 
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The collisional interaction in a Bose condensate represents a non-linearity which in analogy with 
non-linear optics gives rise to unique quantum features. In this paper we apply a Monte Carlo 
method based on the positive P pseudo-probability distribution from quantum optics to analyze 
the efficiency of spin squeezing by collisions in a two-component condensate. The squeezing can be 
controlled by choosing appropiate collision parameters or by manipulating the motional states of 
the two components. 
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I. INTRODUCTION 



In quantum optics it was realized some time ago 
that so-called 'non-classical' states of light, in particular, 
the so-called squeezed states, may outperform classical 
field states in high precision experiments [EI. Moder- 
ately squeezed light has been produced, and several non- 
classical properties have been demonstrated, but classi- 
cal technology has so far left enough space for improve- 
ments that a real practical high precision application of 
squeezed light remains to be seen. Atomic squeezing is 
defined in a way similar to squeezing of light, in the sense 
that the quantum mechanical uncertainty of a physical 
variable, here a component of the collective spin, is re- 
duced by a suitable manipulation of the quantum state 
of the system. 

Due to quite strong interactions between atoms and be- 
tween atoms and light, and due to the long storage and 
interaction times for atoms, the degree of correlations and 
squeezing obtainable in atomic systems exceeds the one 
in light by orders of magnitude. Also, there is a real po- 
tential for practical application of squeezed atoms, since 
current atom interferometers and atomic clocks operate 
at a limit of precision which can only be improved by 
imposing quantum correlations between the atoms. 
Recently there has been a number of proposals for prac- 
tical spin squeezing: absorption of squeezed light 
quantum non-demolition atomic detection pL collisional 
interactions in classical or degenerate gasses TH, [| , photo- 
dissociation of molecular condensates ||, and controlled 
dynamics in quantum computers with ions or atoms [JtJ. 
Since quantum mechanical squeezing implies quantum 
mechanical uncertainties below the level in 'natural' 
states of the system, e.g., the ground state of the sys- 
tems, the analysis of squeezing has to be very precise, 
and in particular one should avoid use of classical ap- 
proximations or assumptions. In the present paper we 
shall investigate the possibilities for squeezing in a Bosc- 
Einstein condensate, using a method which takes the in- 
teractions and the multi-mode character of the problem 



exactly into account. We confirm the validity of results of 
a recent study (^] , and we propose an alternative scheme 
for spin squeezing which relies on a spatial separation of 
atoms in different internal states. 

We consider in this paper a Bose-Einstein condensate 
of two- level atoms in a trap. The dynamics of such a 
system is in the second quantized formalism with creation 
and annihilation operators of atoms in the state i = oor 
b, Tpj , ipi , controlled by the Hamiltonian: 

H = 

i—a,b 

+ g ab ^l(^t(r)Mr)i !a (r)}. (1) 

Here, hi is the single particle Hamiltonian for atoms in 
internal state i and gij is the effective two-body interac- 
tion strength between an atom in state i and one in state 
j. In terms of the corresponding scattering lengths they 
are given by g^ — Airh 2 a,ij /m. 

At temperatures sufficiently below the critical temper- 
ature for Bose-Einstein condensation almost all atoms 
occupy the same single particle wavefunction which to a 
very good approximation can be obtained from a two- 
component Gross-Pitaevskii Equation. There are, how- 
ever, important effects about which the Gross-Pitaevskii 
Equation gives no information: the population statistics 
of the condensate is assumed to be poissonian (or a num- 
ber state ||), it is not treated as a variable which has to 
be determined, and which can be manipulated by phys- 
ical processes. In multi-component condensates, the rel- 
ative populations and coherences of different states are 
important degrees of freedom which require a more elab- 
orate treatment. 



II. TWO-MODE MODEL 

In a first attempt to model the population statistics it 
is convenient for simplicity to assume a separation of the 
spatial and internal degrees of freedom 
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i> a (r) = a(j> a (r) i>b(r) = b(f> b (r). 



(2) 



2 



The separation @ is difficult to justify in general but it 
is certainly reasonable in the initial state that we have 
in mind: A very pure (T=0) single component conden- 
sate is prepared. We then apply a 7r/2-pulse, coherently 
transfcring all atoms to an equal superposition of a and 
b. Then Eq. (||) is fulfilled to a good approximation with 
<i> a = <i>b- 

In the subsequent dynamics, the wavefunctions (f> a (%) 
and 4>b{x) may evolve with time, and we describe this 
evolution with Gross-Pitaevskii Equations 



tonian (W) can be written 



ihdt<t>i = h 



N . |2 
'9n-^m\ + 9% 



N 
~2 



(3) 



where a = b and vice versa. 

The dynamics associated with distribution of atoms 
among the two modes is now studied by rewriting Eq. (Q) 
in terms of the a and b operators: 



^two-mode - 9bb f&TO - yS f S) J dr 3 |0 b | 4 

+ g aa [alataa- y" 1 "^) J dr3 l^a| 4 
+ g ab (WaS - ^b - ^a?a\ J dr 3 |0 h | 2 |0 a | 2 . (4) 

All terms in this Hamiltonian commute with a)d and wb 
and so the interesting dynamics takes place in the co- 
herences between the two internal states. To study these 
coherences it is convenient to define effective internal spin 
operators 



S + — b^a S 1 - = a)b S z = — (i^b — afa 



(5) 



S z represents the difference of total population of states 
a and b, a quantity which can be measured, e.g., by laser 
induced fluoresence in an experiment. S+ and S- or, 
equivalent^, S x = (S+ + SJ)/2 and S y = (S+ - S-)/2i 
represent the coherences between a and b. By suitable 
couplings they can also be turned into population dif- 
ferences and are therefore interesting observables of the 
system. In fact it turns out that the quantity that deter- 
mines the accuracy of a given spectroscopic experiment 
is the ratio of the measured spin-component (the signal) 
to the fluctuations in a component perpendicular to it 
(the noise) ||. For the initial state prepared by a tt/2 
pulse we have N independent atoms all in the same equal 
superposition of a and b. Then the mean spin is in the 
xy-plane and by definition we can take it to be along the 
x-axis. The maximal signal we can obtain is the length 
of the spin S = N/2 and the perpendicular fluctuations 
are then S/2. To improve this "standard quantum limit" 
signal/noise ratio we can introduce correlations among 
the atoms i.e. we can introduce spin squeezing. 

With the internal spin operators of Eq. (^]) the Hamil- 



two- mode 



e(t)N + k(t)S z 
+ E(t)N 2 



D(t)NS z +x(t)S 2 z (6) 



where the time dependent coefficients are given by inte- 
grals involving the time dependent mode functions found 
from Eq.(g). All the ter ms in the Hamiltonian (||) com- 
mute which is an important simplification as the coeffi- 
cients are time-dependent. The term proportional to S z 
will result in a rotation of the spin around the z-axis. 
The N and N 2 terms give different dynamical overall 
phases to states with different total numbers of atoms. 
Such phases are immaterial if we have no external phase- 
standard to compare with and we neglect these terms for 
the purpose of the present work. The A^S'z-term adds to 
the spin rotation with an angle linear in N. It cannot 
be neglected as the direction of the spin (the phase be- 
tween a- and ^-components) can be probed by a second 
tt/2 pulse phase locked to the first tt/2 pulse. The term 
of interest to us is S 2 , the effect of which is well known 
from the work of Kitagawa and Ueda It produces 

spin-squeezing, i.e. it entangles the individual atoms in a 
way that reduces the fluctuations of the total spin in one 
of the directions perpendicular to the average spin. The 
strength parameter of the squeezing operator is given by 

X(t) = Jdr 3 - (g bb |0 fc | 4 +g aa |0 a | 4 - 2g ab |0 a | 2 \c/> b \ 

(7) 

Given the time integral /i — 2 J x(t')dt' of this parameter 
Kitagawa and Ueda provide the analytical expressions for 
the variance of the squeezed spin component 



(AS, 2 

S 
2 



(8) 



where A = 1 — cos 2S_1 /j, and B = 4 sin /i/2 cos 2S_2 fi/2, 
and they specify the direction of the squeezed spin com- 
ponent So = cos QS y + sin 9S Z : 



TT 1 B 

— I — arctan — . 
2 2 A 



(9) 



It is of some interest to use simple analytical approxi- 
mations for x(t) an d we will do so in the specific cases 
studied below. Another approach would of course be to 
obtain x(t) by a numerical solution of Eq. (||). 



III. FULL MULTI-MODE DESCRIPTION 

Equations (§)-(0) are based on a simplifying assump- 
tion. The actual observables of the system are more com- 
plicated to deal with, but we can define a set of operators 



3 



obeying angular momentum commutation relations by 



\ I dr s ^ a + i>Ub) (10) 



1 



J z =\f dr 3 - ftj a ) . (12) 



The total number operator N = J dr 3 ( , 4>'} l '4>a+'4 J l'4>b) com- 
mutes with these three operators and when the two-mode 
approximation applies well, the two-mode and multi- 
mode operators are comparable by the replacement 



J+ = pe iv S+ J_ = pe~ lv S- J z = S z 



(13) 



where pe %v = J dr 3 (j>l(r)(j) a ('r). The factor p takes into 
account that (J x ) and (J y ) vanish unless the atoms in 
state a and b occupy the same region in phase space and 
v is a dynamical phase from the spatial dynamics. 

The two-mode approximation provides an intuitive pic- 
ture of the evolution of the system. It is however not 
easy to justify the factorization (^) and we shall there- 
fore apply an exact method to determine more precisely 
what happens to the mean values and the variances of the 

components of J. The positive P function (P+) [[H] is 
a pseudo-probability distribution giving expectation val- 
ues of normally ordered operator products as c-number 
averages. In our case with two internal states one has 



where 



■■fMt)}-- 



#(*) 

Mt) 



(14) 



and ift 




(15) 



The distribution is not determined uniquely but one par- 
ticular choice obeys a functional Fokker-Planck equation 
which is of course immensely difficult to solve. For nu- 
merical purposes it is much better to translate it to cou- 
pled Langevin equations for the 4 c-number fields ( "wave 
functions"). These equations are "noisy Gross-Pitaevskii 
equations" 

idt 



, , ,, idt /» 
#ip = (-1) [hi + 9% 



(16) 

where a = b and vice versa. In order to treat the inter- 
actions exactly (within the approximation given by the 
form of Eq. ([!])) the noise terms have to be gaussian and 
to fulfill: 



(dW Hl (r,t)) =0 
(dW^(r,t)dW 3 „(r",t')) = 6( 3 \r-7)6(t~t')6^ 



x(-l) 



(17) 



(18) 



On the computer we can simulate the Langevin equa- 
tions to obtain an ensemble of realizations of if). This 
ensemble is a finite sampling of P + and can therefore be 
used to calculate expectation values via the prescription 
( |l4| ) and with a precision limited only by the number of 
realizations in the ensemble. 

Three limitations of the method should be noted: (i) 
The initial state of the system has to be expressed as an 
initial P+. This is trivial for a single coherent state and 
for any mixture of coherent states but it can be compli- 
cated for other initial conditions, (it) The P+ method has 
a notorious divergence problem at large non-linearities. 
This problem sets in after a certain time and is clearly 
noticable in the simulations. We can therefore easily tell 
how long we can trust the results of the method. (Hi) 
Although we have so far written all equations in 3D it 
would be computationally very heavy to simulate a suffi- 
cient number of realizations of Eq. (|16|). In what follows 
we will thus restrict ourselves to ID. It is reasonable to 
assume that this may alter the quantitative results sig- 
nificantly but a ID calculation can be used to investigate 
the validity of the two-mode model which may hereafter 
be applied in 3D with more confidence. 



IV. SPIN SQUEEZING WITH CONTROLLED 
COLLISION STRENGTHS 



Let us first focus on a situation with simple spatial 
dynamics. If we set V a — V& = mui 2 x 2 /2 (ID model) 
and we assume that the values of the collision strengths 
can be controlled so that, e.g., g aa = 56b = ^9ab = 9 the 
spatial dynamics is limited to a slight breathing. The 
spin-dynamics is almost a pure squeezing, that is, the 
mean spin stays in the x-direction. To get an estimate 
of the strength parameter x of Eq. (Q) we find (j> a (t — 
0) = 4>b(t — 0) as the Thomas-Fermi approximation to 
the stationary solution of the GPE with all atoms in the 
a-state. We then have: 



1 



X 



9 



2/3 



2 2 /332/35 N l/3 



m x l 3 u?l 3 . 



(19) 



Choosing g — 0.005fru;ao and N = 2000 we get \ — 
6.1 x 10 _ hw which should give a sizable squeezing within 

a quarter of a trapping periode. ao = \J is the har- 
monic oscillator length in the trap and m is the atomic 
mass. 

In FigJll we show results of both the two- mode approx- 
imation (H) and of the P+ simulation for the parameters 
mentioned above, x was assumed to be constant and of 
the value determined by Eq. (|l^). Jg = cos9J y + sin#J z 
refers to the squeezed component of the spin. The di- 
rection 8 in the yz-plane is determined in the two mode 
model, i.e., from Eq. (0) and it is not independently op- 
timized for the full P+ results. The agreement is seen to 
be surprisingly good considering the crudeness of the es- 
timate of the parameters in the two mode model. Within 
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FIG. 1: Squeezing of the collective spin for favorable in- 
teraction parameters: N=2000 atoms and (g aa , g a b, gbb) = 
(1.0,0.5, 1.0) x 5 x I0~ z hu)a (ID model). + (with errorbars) 
show P+ results, lines show results of the two-mode model. 



FIG. 2: Squeezing due to spatial separation of the two inter- 
nal state potentials by a;o = 3ao. The total number of atoms 
is 2000 and the interaction strengths are (g aa , gab, gbb) = 
(1.03, 1,0.97) x 5 x 10- 3 hua (ID model). 



the uncertainty of the positive P results (the (AJ|) data 
are presented with errorbars in Fig. |l|) the noise sup- 
pression is seen to be almost perfect. These calculations 
thus confirm the results of S0rensen et al. || where the 
simple two-mode approach was supplemented by another 
approximate method. 



SPIN SQUEEZING WITH CONTROLLED 
MODE FUNCTIONS OVERLAPS 



In the experiments on the \F = 1,771/ = —1) and 
\F = 2, nif = 1) states in 87 Rb |13j the scattering lengths 
and thus the interaction strengths are actually in propor- 
tion g aa : g a h : g^h — 1-03 : 1 : 0.97. This is far from ideal 
conditions for squeezing as we can see from Eq. (Q) : \ = 
when in addition cj> a = <ph as is the case initially. To pro- 
duce a sizable squeezing effect we propose to make the 
two mode functions differ (see also |u| ) . This is achieved 
by applying different potentials to the two internal states 
which has actually already been done for magnetically 
trapped Rb making use of gravity and different magnetic 
moments of the two internal states |0|. If V& is dis- 
placed from V a the fr-component created by the initial 
7r/2-pulsc will move away from the a component thereby 
reducing the overlap Jdx\4> a \ \4>b\ and increasing \- It 
is remarkable that the squeezing then takes place while 
the two components are away from each other and are 
therefore not interacting. 

To demonstrate the accomplishments of the scheme de- 
scribed above we have chosen simply to displace Vh by a 



certain amount xq from V a . In this model both potentials 
are still harmonic and of the same strength. The spa- 
tial dynamics is now more complicated but it can be ap- 
proximated by the solutions to coupled Gross-Pitaevskii 
equations. To get a rough idea of the evolution we can 
use the well know evolution of a displaced ground state 
wavefunction in a harmonic trap. We then get 



iu — Xq(\ — cos(wt))/2 ~ix Q sin(aji)/2 



and 



pe 



\^>b\ 2 \^>a\ 2 dx 



-Xo(l-cos(wt))/4 



(20) 



(21) 



If we choose (g aa , g a b, gbb) = (1-03, 1, 0.97) x 5x 10~ 3 hwao 
and xq — 3ao this model gives an order of magni- 
tude estimate for the integrated strength parameter of 
J 21 '" x{t)dt = lQ~ 2 fiw when the two components are 
again overlapped. For 2000 atoms this corresponds to a 
reduction of the uncertainty in the squeezed spin compo- 
nent by roughly a factor of 10 according to Eq. ||. 

In Fig]| we use these estimates to analyze the P+ sim- 
ulations, i.e., we plot the expectation value of the spin 
component predicted to be maximal, ( J v ) , and the vari- 
ance of the perpendicular component predicted to be 
squeezed, (AJg). As can be seen, after a fast drop a 
large fraction of the original mean spin is recovered in 
( J„) when the two wave packets are again overlapped. 
This confirms our prediction of v and it implies that a 
sizable signal can be obtained in an experiment. At the 
same time the variance in the predicted perpendicular 



5 



component is strongly suppressed confirming our predic- 
tion for 9 and implying that the noise in the experiment 
can be significantly reduced below the standard quantum 
limit. 



VI. DISCUSSION 

We have presented the first exact numerical studies of 
spin squeezing in two-component condensates. As in our 
previous work on squezing of a single component atomic 
field operator [|| we have applied the positive P distri- 
bution which is well suited for studies of transient, short 
time behavior of interacting many-body systems. This 
formulation involves simulations, and due to the formal 
separation of c-number representations of creation and 
annihilation operators, 'non-physical' results may occur 
such as negative/complex atomic densities, and few in- 
stants with negative variances as depicted in Figs, [j] and 
||. For long times we cannot exclude that slight impreci- 
sion due to discretization of space and time contributes 
to these results. The Monte Carlo nature of the method 
makes it somewhat tedious to check rigorously for dis- 
cretization errors but for resonable numbers of realiza- 
tions in the ensemble sampling errors will dominate any- 
how. 

The conclusion of this paper is that the predictions of 
the simple two-mode model reproduce the full multimodc 
result also quantitatively. This means that the spatio- 



temporal dynamics and the population dynamics couple 
the way they should to produce spin squeezing but the 
resulting entanglement between spatial and internal de- 
grees of freedom is small enough that purely internal state 
observables show strong squeezing and multi-particle en- 
tanglement. In particular the calculations of Fig. || re- 
veal that quite significant distortions of the distributions 
when the components separate and merge do not prevent 
sizable spin squeezing. 

Although we claim quantitative agreement above it 
should of course be realized that this is within the ID 
model. To use the two-mode model to describe a real 
experiment the spatial dynamics which enters the inter- 
nal dynamics via Eqs. ([?]) and ( |l3| ) must be treated with 
some accuracy. Apart from the 3D aspects it should be 
noted that e.g. in some experiments on two-component 
condensates [[L4[ the displacement of the potentials are 
accompanied (and in fact due to) different strenghts of 
the potentials. This adds to the complexity of the dy- 
namics and will of course be important if the aim is pre- 
cise quantitative predictions and not a proof-of-principle 
analysis as offered in this paper. 

More tests need to be carried out for other kinds of 
processes, but the present study suggests that the sim- 
ple two-mode description provides good predictions for 
the many-body dynamics of spin squeezing. Hence other 
ideas, e.g., for reducing the fluctuations in the total num- 
ber of atoms in a condensate or in out-coupled atom laser 
beams, may be reliably based on the terms of Eq. (0). 
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